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We extend the Weak Crystallization theory to the case of metallic alloys. The additional ingredient 
- itinerant electrons - generates nontrivial dependence of free energy on the angles between ordering 
wave vectors of ionic density. That leads to stabilization of FCC, Rhombohedral, and icosahedral 
quasicrystalline (iQC) phases, which are absent in the generic theory with only local interactions. 
The condition for stability of iQC that we find, is consistent with the Hume-Rothery rules known 
empirically for majority of stable iQC; namely, the length of the primary Bragg peak wavevector is 
approximately equal to the diameter of the Fermi surface. 


Introduction. Crystallization is probably the most fa¬ 
miliar but one of the hardest to analyze phase transi¬ 
tions. The workhorse of the theory of phase transitions, 
the Ginzburg-Landau theory |T] , cannot be easily applied 
to many of the crystallization transitions since they tend 
to be strongly first order, i.e., the order parameter expe¬ 
riences a large jump at the transition. 

Crystals are best characterized in reciprocal space, 
where the onset of long-range order is signaled by the 
appearance of resolution-limited Bragg peaks. The in¬ 
tensity of the Bragg peaks reflects the density distribu¬ 
tion in a material - for smooth density modulations, as 
in the case of liquid crystals, only few harmonics of prin¬ 
cipal peaks located on a momentum shell of radius qo 
are needed to fully describe the state. Then, the inten¬ 
sity of principal harmonics is the order parameter, and 
when it is small near the transition (relative to the av¬ 
erage density), the application of the Ginzburg-Landau 
theory is justified. In atomic crystals, however, density 
is highly peaked at the equilibrium positions of atoms, 
and the number of relevant Bragg peak harmonics scales 
in proportion to the ratio of the unit cell size to the 
atomic size (smeared by thermal and quantum fluctu¬ 
ations). In a typical crystal, the thermal fluctuations of 
atoms are 15 — 30% of the lattice spacing at the melting 
transition |2]; therefore, to accuratley describe the tran¬ 
sition, multiple harmonics of qo are required. The ap¬ 
pearance of strong modulation immedaitely at the phase 
transition, with multiple Bragg peaks forming recipro¬ 
cal lattice, is the signature of strongly first order transi¬ 
tion. A special case of a crystalline solid is a quasicrystal , 
where atoms lack simple spatial periodicity; yet, in the 
reciprocal space, resolution-limited Bragg peaks appear 
in a self-similar arrangement inconsistent with crystallo- 
graphically allowed point-group symmetries [3Hz] . 

Weak Crystallization theory [8] applies Ginzburg- 
Landau machinery to crystallization by assuming that 
only the Bragg peaks on a single momentum shell are 
significant enough to affect energy. While most natu¬ 
rally applicable to liquid crystals, the theory has been 
used to predict ubiquity of Body Centered Cubic (BCC) 
crystals near crystallization temperature [9], to study ef¬ 
fects of fluctuations m, and, with some modifications, 
to address the problem of stability of quasicrystals m 


ITS] . As such, it has been a useful symmetry-based tool 
to study the crystallization transition, even beyond its 
immediate range of validity. However, in its standard 
spatially-local form, weak crystallization theory is inca¬ 
pable of obtaining many of the experimentally relevant 
crystalline states, such as simple cubic, rhombohedral, or 
Face Centered Cubic (FCC), and the heuristic modifica¬ 
tions of it have not been microscopically justified. 

Crystals are often (qualitatively) separated into classes 
based on the dominant type of interaction that holds 
them together: ionic, covalent, molecular, and metallic. 
The crystal structure depends on a variety of details, such 
as ionic charge and electronic orbital structure, etc. It 
is therefore remarkable, that in the case of metallic crys¬ 
talline alloys simple empirical rules exist. These rules 
were identified by Hume-Rothery m who has found that 
metallic alloys are particularly stable when in addition to 
the requirement that atoms be of similar size and elec¬ 
tronegativity, the value of the average valence per atom 
(“e/a” ratio) be close to certain “magic” values, which 
depend on the crystal structure. The optimal e/a ratio 
has been argued to be associated with a particular geo¬ 
metrical matching condition, when the itinerant (nearly- 
free) electron Fermi surface “just crosses” the boundary 
of the first Brillouin zone m- Regardless of interpre¬ 
tation, this observation points to an important, if not 
decisive, role that itinerant electrons play in determining 
the crystal structure. This is indeed not surprising given 
that electrons are an effective mediator of long-ranged 
and multi-body interionic interactions. 

Significantly, majority of stable quasicrystals have 
turned out to be Hume-Rothery alloys [T8], i.e., they are 
stable for narrow ranges of e/a. Despite nominally large 
conduction electron concentration, their electrical and 
thermal conductivities are exceptionally low PI , con¬ 
sistent with strong scattering around the Fermi surface. 
In analogy to regular crystals, attempts have been made 
to construct a theory of quasicrystals accounting for the 
Hume-Rothery rules by perturbatively including electron 
scattering on quasiperiodic ionic potential [20] . Just as 
in the case of regular crystals, such approach is problem¬ 
atic (see Appendix 0. There have also been attempts 
to understand the formation of quasicrystals in terms of 
Weak crystallization theory P3HI5]; however, there have 
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been no microscopic justification of the theoretical as¬ 
sumptions, and in particular these theories cannot ac¬ 
count for the fact that most of quasicrystals follow the 
Hume-Rothery rules. 


Here we extend the Weak Crystallization theory to 
metallic systems. We do so by explicitly introducing itin¬ 
erant electrons that couple to the ionic density, and in¬ 
tegrating them out. We find that interionic interactions 
generated by electrons qualitatively modify the generic 
weak crystallization theory, stabilizing FCC, rhombohe- 
dral, and, notably, icosahedral quasicrystal (iQC) states. 
The Hume-Rothery rules emerge from the interplay of 
two length scales - the preferred interionic distance, 1/go, 
and the Fermi wavelength of itinerant electrons, 1 jkp. In 
particular, for iQC we find go ~ 2 kp, consistent with em¬ 
pirical observations m- 


In our approach we explicitly calculate the electronic 
contributions to the quadratic, cubic, and quartic in den¬ 
sity terms in the Ginzburg-Landau (GL) energy. In this 
regard it is analogous the direct derivation of Free energy 
in the cases of superconductivity m and charge density 
waves [22]. Of crucial importance is the strong (singular 
in the limit of zero temperature) dependence of quartic 
term on the angle between the Bragg wave vectors. At 
some angles, which can be tuned by changing the size of 
the electronic Fermi surface, the repulsion between differ¬ 
ent Fourier components of ionic density can be reduced 
or even turned into attraction. For a generic Fermi mo¬ 
mentum kp, this effect favors formation of rhombohedral 
states, with three ordering wavevectors with identical an¬ 
gles a m in(kF ) between them. For more finely tuned kp, 
the minimum interaction can be reached at angles a m i n 
(and equivalent n — <a m ^ n angles) that favor simple cubic, 
FCC, or iQC. It is interesting to note that FCC state 
additionally benefits from a noncoplanar fourth order in¬ 
teraction terms that can be always made attractive by 
appropriate choice of the order parameter signs, expand¬ 
ing the domain of its stability. 


Weak crystallization theory and its extension to met¬ 
als. In what follows, we shall keep only momenta of 
length g 0 ; he., we shall make the ansatz p(x) = p 0 + 
X!|k|=qo R e [PkC k x ] for ionic density. As discussed above, 
this ansatz, which is central to “weak crystallization” 
theory [8], is strictly valid only where the crystalliza¬ 
tion transition is weakly first-order and its latent heat 
is small. Outside this regime, our results will not be 
quantitatively accurate; nevertheless, we expect them to 
provide guidance as to what kinds of crystal structures 
are favored. 


We proceed by writing down a general Ginzburg- 
Landau (GL) free energy functional, T = Fq + Hq + V, 
where 




FIG. 1. (a) First shell of reciprocal lattice vectors for the case 
of FCC real space (BCC reciprocal) lattice vectors. Four vec¬ 
tors point from the center of a tetrahedron to its vertices, (b) 
Example of a scattering path between electronic momentum 
states induced by the ionic densities modulated at wavevec¬ 
tors in (a). This process makes contribution to ie({qi}) in 
Eq. Q. (c) Feynman “box” diagram for the fourth order 
contribution to the GL free energy. Lines correspond to elec¬ 
tronic Green functions; vertices to ionic densities with a given 
wavevector. (d) An example of coplanar contribution to free 
energy that only contains two pairs of ±q^ [a contribution to 
u(qo,aij)\. As shown, corresponds to the “resonance” con¬ 
ditions satisfied: tree or more electronic momenta are on a 
great circe of the Fermi surface. 


p o = Yi r 0?)M 2 + ^7 Yi PgiPtoPg^C^Zvi) 

U O' Q* 

PqiPqiPqzPqi^C^jQi) ( 1 ) 

ffo = E k [£(fc)-M]Tk (2) 

V = ^2 kcl V ^Pl C k C k-q ( 3 ) 

Here, Fq describes the physics of ions and core elec¬ 

trons in the absence of itinerant electrons. The minimal 
(“local”) assumption that is commonly made is that in¬ 
teractions A3 and A4 are simply constants. However, as 
interatomic interactions set a preferred length-scale for 
crystallization even in the absence of conduction elec¬ 
trons, this length-scale is introduced into the second- 
order term, via the weak-crystallization form r(g) = 
r o+x(M — go) 2 - As already stated, we will restrict our at¬ 
tention to density modes that are precisely at q = go [23]. 
The second term, iLo, describes the itinerant electrons: 
for simplicity we shall treat these as noninteracting. The 
third term, V, describes the interaction between itinerant 
electrons and atoms. As we are only concerned with den¬ 
sity modulations satisfying g = go, and the interaction is 
assumed to be spherically symmetric, we can parameter¬ 
ize the interaction strength entirely by its Fourier compo- 
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FIG. 2. Electronic contribution u(a) to 4th order term in 
Ginzburg-Landau energy functional Eq. Q. Temperature 
is T = 0.1 Ef- Black lines mark location of zero-temperature 
singularity, Q/kF — 2 cos f. 


nent at momentum transfer q 0 , viz. v = v(qo). Thus we 
need not make any assumptions about screening of the 
Coulomb interaction. The kinematic constraints ^ ., in 
combination with the q = qo restriction, strongly limit 
the number of allowed terms. Namely, the cubic term is 
only non-zero for triplets of forming equilateral trian¬ 
gles, and thus favors hexagonal and BCC crystal struc¬ 


tures [ 9 ]. The quartic term obtains generically from com¬ 
bining ±q i with =tqj. It can also appear in the situation 
when four q^ form a non-copanar quadrilateral [e.g., the 
geometry in Fig. |TJa)] . 

Electronic contribution to Weak crystallization energy 
functional. We now integrate out the conduction elec¬ 
trons to arrive at a description that is purely in terms of 
the atomic densities. The procedure is analogous to the 
derivation of Ginzburg-Landau functional for supercon¬ 
ductivity or a charge density wave states. The difference 
is that the ionic density order parameter, a priori, can 
have an arbitrary number of components, and the en¬ 
ergy functional should be able to predict not only the 
magnitude of the order parameter, but also the number 
and orientation of its components. The latter determines 
the type of crystalline state. 

As the free energy functional T = Fo + Ho + V is 
quadratic in fermion operators, we can integrate out the 
fermions; this allows us to write the partition function 
purely in terms of ionic densities, as Z = exp[— / 3 (Fq — 
AT 1 )] where Fq is defined in Eq. 0 and AF is given 
by the following perturbation series, which we have re¬ 
summed using the linked cluster theorem: 


a f = 


p 


/ 


dT 1 ...dT n (T T V(T 1 )...V(T n )) c 


Explicit expressions are given in Appendix [B] We now 
expand A F to quartic order in the bosonic densities; this 
yields, for the free energy functional F = Fq + AF: 


F — IPql ^ 3 (^ 0 ) PqiPq2Pq3^(^~^ Qi) 

+ \ H q ^ qj J A4 + “(“«)] K| 2 K I 2 + \ + W(0)]|p q j 4 + [A4 + re({qi})]p qi Pq 2 Pq 3 /)q 4 ( 4 ) 


The symbols and indicate summation over 
unique triangles and non-planar quadrilaterals of qi; otij 
is the angle between vectors q^ and q^. 

Numerical results. Figure [2] shows u(a) for vari¬ 
ous values of qo at T = O.li^i?. Already at this not 
very low temperature certain features become apparent. 
For qo/kp ~ a minimum in u(a) develops around 
a = 7 r/ 2 , which then splits into two minima for larger 
values of q^. In the limit of zero temperature, a sin¬ 
gularity develops along the line qo/kp = 2cos^. Ge¬ 
ometrically this condition corresponds to the configura¬ 
tion when three momenta connected by scattering off the 
ionic order parameter, k, k + q x and k + q2 can all si¬ 
multaneously be on a great circle of the Fermi surface 
(Fig. [IJi). Near this line, the vertex is repulsively diver¬ 
gent for smaller qo and attractively divergent for larger 
qo as T —»• 0 . This singular behavior is a four parti¬ 
cle/hole analog of the particle-hole divergence in ID that 


drives Peierls instability. Naturally, such a strong angu¬ 
lar dependence of u(a) at temperatures much lower than 
Ep can influence the energetic balance between different 
crystalline phases. It should be noted, that the angular 
dependence is a result of sharply defined Fermi surface; 
and temperatures comparable or higher than the Fermi 
energy it becomes smeared out. 

The non-coplanar terms w({qi}) are less generic than 
u(a) since they require four distinct wavevectors to add 
up to zero. A case where these terms are important is the 
FCC crystal, whose first Bragg shell (the set of shortest 
symmetry-related reciprocal lattice vectors) is comprised 
of eight vertices of a cube; hence there are two non-trivial 
quadruplets of wavevectors that correspond to the ver¬ 
tices of two tetrahedra (see Fig. [l]). The significance of 
non-coplanar terms is that they can always be made to 
lower energy by appropriate choice of the relative signs 
of constituent p qi . For the momentum-independent in- 
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FIG. 3. Temperature and qo dependence of noncoplanar elec¬ 
tronic contribution w(qo) to 4th order term in Ginzburg- 
Landau energy for FCC crystal. Black vertical line marks 
location of zero-temperature singularity, qo/kF = 2cosf, 
where a = arccos(l/3). 


teractions this does not change the fact that stripe state 
has the lowest energy within Weak Crystallization theory. 
However, inclusion of the electron-induced interaction 
can change he situation dramatically. In Figure [3] we plot 
re({qi}) as a function of qo for FCC. We find that it has 
features similar to u(a)\ namely, when qo/kp = 2 cos 
with a t = arccos(l/ 3 ) the tetrahedral angle, w diverges 
as T — y 0 . This enhanced interaction is the cause of a 
large region of stability of FCC phase that we find. 

Phase diagram. To construct the phase diagram, we 
first consider the set of variational states that contain N 
pairs of =Lqwhere all q^’s are symmetry related, and 
hence have exactly the same set of neighbors. Then, all 
the Fourier amplitudes are identical, |p q . | = p and the 
free energy is 


F = Nr 0 \p\ 2 +^ 


^u(aoj) + 
j/o 


2 Mq 
N 




M 4 , 


where we redefined u = u + A 4 and w = w + A 4 for 
compactness. We assumed that all M\j quadruplets have 
the same re({qi}) (the case for FCC) and that vectors 
qi do not form equilateral triangles, and hence the cu¬ 
bic invariant that could stabilize BCC (FCC reciprocal) 
crystal does not contribute (the latter assumption should 
become valid for sufficiently large negative 7*0). Now it 
only remains to minimize the energy to obtain, 


H 2 


r o 


and 

F = _ - _ (5) 

It is important to note that the pure electronic vertex 
u is negative (attractive) in a wide range of qo and <a, 


which taken by itself could cause an absolute instability. 
In this regime, one cannot truncate F at fourth order, but 
must include higher-order terms in the GL expansion to 
find stable equilibrium states. However, the structureless 
local interaction A4 restores stability while maintaining 
the strong angle-dependence of the interactions. 

The results of our analysis are presented in Figure [ 4 ] 
We find only four stable phases: rhombohedral, striped, 
FCC, and iQC (i.e., icosahedral quasicrystal). The 
other symmetric variational states we explored are al¬ 
ways higher in free energy than these (see lower panel 
of Figure [ 4 ] for energy comparison and Appendix [C] for 
details of the variational states). The overall shape of 
the phase diagram can be understood as follows. When 
the structureless interaction A4 is absent or too weak, 
as noted above, the free energy can become unstable 
at quartic order. On the other hand, when A4 is dom¬ 
inant, the electron-induced interaction can be ignored, 
and we recover the standard weak-crystallization result 
that the equilibrium state is striped (or smectic). When 
we are far from the matching condition qo ~ 2 kp, or the 
temperature is relatively high, the interactions are not 
strongly angle-dependent, and these are the two domi¬ 
nant possibilites. On the other hand, when the struc¬ 
tureless and electronic contributions are of similar mag¬ 
nitude and the temperature is low, the angle-dependence 
of the electron-mediated interaction stabilizes nontrivial 
crystalline phases. 

The most significant qualitative feature of the phase 
diagram for intermediate values of A4 is the dominance 
of FCC and rhombohedral phases, with balance shift¬ 
ing in favor of FCC at lower temperatures. The rea¬ 
son for this trend is that FCC has two appealing fea¬ 
tures - (1) it has only one inter-q angle (up to 
7r — dij ), and ( 2 ) it has lower energy due to the pres¬ 
ence of non-coplanar 4 th order terms. Rhombohedral 
crystal only has former feature, and thus only becomes 
competitive when the angle <a m * n that minimizes u(a) is 
sufficiently far from the tetrahedral angle. Surrounded 
by the FCC and rhombohedral phases is the iQC phase. 
The key advantage of iQC phase is that it has large num¬ 
ber (six) of ±q i pairs, all separated by the same angle 
OL{ = 2sin -1 (7 2 + l) -1 ~ 63 . 43 ° (7 is the Golden mean). 
Even thought iQC can not benefit from the non-coplanar 
energy terms, when the optimal angle <a m * n is close to the 
ai, iQC can beat both FCC and Rhombohedral. Finally, 
for large A4 we recover the Stripe phase predicted by the 
original featureless Weak crystallization theory. 

Apart from the FCC phase with its non-coplanar terms 
in energy, the phase diagram can be understood as fol¬ 
lows. For the states with only one non-trivial inter-q an¬ 
gle a m £ n , the denominator in Eq. © is ii(c^mm) T [^(0) 
2u(a m i n )] / N . Thus, if the second term is positive, it fa¬ 
vors large N ; it it is negative, then N = 1 . Note that 
for 7/(0) = 2 u(a) states with all possible TV’s are energet¬ 
ically degenerate (can be seen in triple crossing point in 
Figure flower) at qo/kp ~ 2 . 05 ). 

In construction of the phase diagram we have only 
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FIG. 4. Variational phase diagram for T = 0 . 2 Ef (top) and 
T = 0.1 Ef (middle) as a function of local repulsion strength 
A 4 and the ionic ordering vector g 0 . At small values of A 4 
(region marked Unst), the fourth order terms in GL for one 
or more of the variational state becomes negative, signaling 
the need to consider higher order stabilizing terms. Lower 
panel shows Energy comparison of different variational states 
for a fixed value of A 4 at T = 0.1 Ef- 


considered highly symmetric states. We now discuss 
possible deviations from these assumptions. First, we 
can ask whether the highly symmetric crystal states are 
stable with respect to “Bragg-fractionalization,” namely 
whether it may be beneficial to split Bragg peaks into 
multiple nearby ones. From the fact that u(a 0) = 
ffo, which can be explicitly demonstrated for electron- 
mediated and local interaction, for uq > 0, lack of frac- 
tionalization follows trivially (see AppendixjD]). The next 
possibility is a distortion of peaks from symmetric posi¬ 
tions. Clearly this is not a concern for Rhombohedral 
state, but could be for iQC and FCC. Here we specifi¬ 
cally ask whether iQC will remain stable even if a m i n is 


not exactly (see Appendix 0. Due to its high sym¬ 
metry iQC cannot naturally distort, unlike, e.g., Rhom¬ 
bohedral state. To answer this question, we have ex¬ 
panded the interaction energy around the symmetric iQC 
state. We have found that if v!(on) < —(2/3 )u"(ai), 
then iQC spontaneously distorts into a lower symme¬ 
try state, i.e. a distortion could occur if <a m ^ n > cq 
(“compresses springs”). This criterion also shows that if 
u(a) is sufficiently smooth, as it is for temperatures not 
very much smaller than the Fermi energy, then undis¬ 
torted iQC should in fact be quite stable. Indeed, ex¬ 
panding around A, we find that criterion for instability 
is OLmin ~ > 2/3, i.e., the minimum is at least 40° 

away (above) from the icosahedral angle. This estimate 
is based on the assumption of smoothness of u(a), which 
is violated at temperatures much below the Fermi tem¬ 
perature. Thus, for quasicrystals that form under such 
conditions, there is a possibility of distorted iQC, as well 
as a structural transition from perfect to distorted iQC 
as a function of temperature. 

We have also explored possible ordered states in an 
alternative fashion by applying unconstrained stochastic 
minimization of the free energy functional. To simplify 
simulations, we neglected the cubic and non-coplanar 
quartic terms and thus cannot fully capture FCC and 
BCC phases; however, the advantage of this method 
is that it provides an unbiased treatment for arbitrary 
multi-g states that are not required to possess any spe¬ 
cial symmetries. We start from random configuration of 
several hundred components p qi with on a sphere or 
radius go. We then iteratively minimize energy by ran¬ 
domly selecting p q . and changing its value and position 
on go sphere in the direction of decreasing energy. The 
minimization results are consistent with the variational 
phase diagram in Figure [4] (modulo underestimating the 
stability of FCC). Due to the stochastic nature of the al¬ 
gorithm, however, it sometimes converges to other states. 
In particular, in the region of stability of iQC, the final 
state is rather commonly the decagonal state, which is 
approximately the iQC state with one pair ±g^ removed. 
This state is and example of a 2d quasicrystal - it is 
periodic along one axis and quasiperiodic in the plane 
perpendicular to it. Even though the energy of this state 
is very close to the iQC, we have not observed it ever 
to be lower in energy than the perfect iQC (consistently 
with Figure [4] (lower)). The energy difference is never¬ 
theless sufficiently delicate, so one cannot rule out that 
for modified conditions decagonal state may appear as 
the lowest energy state in the phase diagram. 

The conjecture that stability of 3D quasicrystals is as¬ 
sociated with “bond-orientational order” that favors spe¬ 
cific inter-q^ angles within Weak Crystallization theory 
has been previously expressed by Mermin and Troian [T2] 
and Jaric [13]. In Ref. [12] an auxiliary field was intro¬ 
duced to generate preferred inter-q^ angle, however, no 
physical justification was given as to the nature of this 
field. The key result of our work is that itinerant elec¬ 
trons play the role similar to the auxiliary field postu- 
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lated in [12]. On the experimental side, it has been found 
that the optimal e/a ratio observed in quasicrystals cor¬ 
responds to the approximate matching between the qua- 
sisrystalline gtmsz-Brillouin zone and the electronic Fermi 
surface, that is, the length of the dominant Bragg wave 
vector equals the diameter of the Fermi surface, 2 kp- 
This is indeed what we find. 

In conclusion, we have analyzed the effects of electron- 
ion interactions on crystallization transition within Weak 
Crystallization theory. We found that the angular- 
depended multi-ion interactions induced by electrons can 
lead to stabilization of such empirically common but elu¬ 
sive, within the standard theory, states as Rhombohe- 
dral, FCC, and icosahedral quaiscrystals. The stability 
conditions can be recast in terms of Hume-Rothery rules 
connecting primary ionic ordering wave-vectors and the 


size of the electronic Fermi surface. Our results are ob¬ 
tained within the assumption that the cubic invariants 
are less relevant than the quartic ones, i.e., at temper¬ 
atures sufficiently lower than the temperature of mean 
field transition (ro = 0). Near the transition, more care¬ 
ful analysis of fluctuations is required m, which will be 
the subject of future work. 
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Appendix A: Other approaches to crystallization 
theory and their limitations. 

The most common approach to determine the lowest- 
energy crystal structure is based on variants of micro¬ 
scopic density functional theory, which specifies atoms 
with their electronic shells and attempts to optimize their 
spatial arrangement [24] . Due to computational complex¬ 
ity, this is a variational approach that can effectively treat 
only periodic arrangements of atoms. Near the melt¬ 
ing transition, application of this method becomes diffi¬ 
cult since atoms in a liquid lack spatial periodicity. In 
that regime, methods combining density functional the¬ 
ory with molecular dynamics are applied, but only with 
limited success [25] • 

Periodic approximants to quasicrystals have also been 
studied by density functional theory [26] : application of 
this method, however, requires a very large number of 
atoms to be explicitly considered and optimized for the 
approximants energies to provide a good estimate for 
quasicrystals, even away from the melting transition. 

Another, semi-microscopic, approach is based on the 
Peierls instability-type arguments. There, one studies 
the features in the electronic susceptibility and attempts 
to use its anomalies as a predictor of stable phases. This 
approach is problematic in the case of 3D alloys, as can 
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be easily seen. We would like to have an unbiased predic¬ 
tor of an ordered state; therefore, the only starting point 
possible is free electron Fermi sea coupled to featureless 
(constant) ionic density. In ID, electronic susceptibility 
diverges at 2k p at T = 0, which leads to a density in¬ 
stability at this wavevector; this is the origin of charge 
density waves in many quasi-ID materials. In contrast, 
in 3D, free electron susceptibility is maximized at zero 
momentum and at 2k p only has infinite first derivative 
(the cause of Friedel oscillations of electron mediated in¬ 
teraction [27]). However, this is insufficient to cause in¬ 
stability in ionic density - theory would predict that the 
instability should occur at zero wavevector, i.e., at uni¬ 
form density. Moreover, the (quadratic) term in the GL 
theory that is proportional to the electronic susceptibil¬ 
ity only includes a single density modulation, and thus 
cannot discern between orderings that contain multiple 
wave vectors. 

A way to go beyond the quadratic energy approxima¬ 
tion is to include the ionic modulation non-perturbatively 


in electron dispersion m It has been argued this way 
that for a given crystal structure, the electronic energy 
is minimized when the Fermi surface “just crosses” the 
Brillouin zone boundary. This naturally corresponds to 
crystal-specific optimal e/a ratios, and thus appears to be 
consistent with the empirical Hume-Rothery rules. Ap¬ 
plication of this approach to discriminate between ener¬ 
gies of different crystalline and quasicrystalline states is, 
however, problematic, as it presupposes the knowledge of 
the amplitude of the periodic lattice (pseudo-)potential, 
which is different for different crystals. Since the energies 
of various states are typically rather similar, the uncer¬ 
tainty in the potential makes such approach unreliable. 

Appendix B: Details of electronic corrections to Free 
energy 

Integration of electronic degrees of freedom leads to the 
following corrections to the ionic free energy functional: 


A i J dridT2(T T E V qi p qi c\ i+qi cki I T1 v q2 p q2 c\ 2+q2 ck 2 1 T2 ) conn (Bl) 

ki,qi,k 2 ,q2 

= ^0- J dT 1 dT 2 G p (T 2 - ri)Gp_ g (n - r 2 ) = 00- E G p(. u ’n)Gp- q (ui n ), (B2) 

AF (3) = S(£ qi ) £ GpM G p _ qi (io n ) Gp _ qi _ q2 (to n ), (B3) 


^^(4) _ V qo kq l Pq 2 Pq3 Pq± ^ ^ 

^ uj n ,k 


Here, G p (u n ) = (iu n - e p ) 1 . 

As one can see, the new terms in the Ginzburg-Landau 
functional have the form similar to those already con¬ 
tained in F 0 (Eq. ([!])). The second order term A 
only serves to redefine go and hence will be of no inter¬ 
est to us. The prefactor of the cubic term becomes a 
function of go- We find by numerical integration that the 
electronic contribution to this term is non-singular in the 
limit of zero temperature, and thus it does not introduce 
any qualitatively new features relative to those already 
in Fq. 

The 4th order correction AF( 4 ) can diverge if certain 
geometric conditions are satisfied (see Figure 1 of the 
main text). Thus we concentrate here on this term only. 

The first type of 4th order term that we will consider is 
self-interaction. Self-interaction is generated by box dia¬ 
grams with momentum transfers (qi, qi, —qi, —qi) (type 
A) and (qi, —qi, qi, —qi) (type B). The combinatorial 
multiplicities of these diagrams can be calculated as fol¬ 
lows. At every vertex of the box diagram we can place 
p± qi. For A type, there the sign pattern has to be such 


)G p - qi (iJ n )G p - qi - q r 2 (uJ n )G p - qi - q2 -q 3 (id n ). (B4) 


that same signs are adjacent, while for B - interlaced. 
There are 4 ways to place two adjacent ++ in 4 boxes. 

(+H-), (—bH—), (-1—b), (H-b)- Hence A type 

has multiplicity 4. For B type, there are only two distinct 

ways to arrange: (H-1—), (—}-b), and multiplicity 

is 2. Therefore, the self-interaction goes as 

J2(4A + 2B)lpJ 4 . 

There are two types of mutual interaction diagrams: 
those that contain only two pairs of ±q^ (“coplanar” 
diagrams), and those that contain four distinct q^’s 
(“non-coplanar” diagrams). Coplanar diagrams depend 
only on one angle between qi and q2 (for ft = 0 we 
get self-interaction). There are three distinct contribu¬ 
tions to AF^ 4 ), which come from the following arrange¬ 
ments or momenta around the box diagram: AF^ : 
(qi,-qi,q 2 ,-q 2 ) (type Vi), A Fp : (qi,-qi,-q 2 , q 2 ) 
(type V 2 ), and AF 3 (4) : (qi, q 2 , -qi, -q 2 ) (type D). 
Their combinatorial multiplicities are as follows: 








Vi + V 2 - 4 ways to place qi, 2 ways to place — qi next 
to it (PBC), 2 way to place ±q 2 : total 16. Hence there 
are 8 diagrams of each type. 

D: 4 ways to place qi, 2 ways to place ±q 2 - Total 
multiplicity is 8 . 

Therefore, the mutual interaction term is 


EioW + 8 F 2 + 8 J D)| / 0 qi | * 2 K i | 2 (B5) 

= E^(4U +4V 2 + 4D)K i |> q y. (B 6 ) 


Notice, that in the limit q^ —>> q j, V\ —>> F>, and 
(V 2 ,D) —» A Hence, going back to the original nota¬ 
tion in terms of u(a) , we find that the full 4th order GL 
term is 

sf(4) = \^2 u { a ij)\p^\p^\ 2 + \Yj u °^ 2 

i^j i 

where u(a) = 8 (Vi + V 2 + D) and uq = u(a = 0) (it can 
be shown explicitly that the limit a 0 is continuous at 
finite temperature). 

The interaction functions can be obtained by a mixture 
of analytical and numerical integration. The frequency 
summations could be performed with the help of contour 
integration, 


ry y 1 — V v° l^i l 2 l^ 92 1 2 n ^( e 1) , tif{c 4) 

1 ’ 4 ^ ( Cl - C2 ) 2 (ci - C 4 ) ( C 4 -C 1 )(C 4 -C2 ) 2 

_ ripje 2) _ np(^ 2 ) _^_ ^j?( e 2) 

(e2 — Ci) 2 (e2 — 64) (e2 — ei)(e2 — C4) 2 (e2 — Ci)(e2 — 64) ’ 

and 

jj _ v to\P qi 1 2 1 ^ 21 2 ^ 2 ) 

4 “ (ei - e 2 )(ei - e 3 )(ei - e 4 ) (e 2 - ei)(e 2 - e 3 )(e 2 - e 4 ) 

_ ^(^ 3 ) _^_ ^f(€4) _ 

(^3 - ei)(e 3 - e 2 )(e 3 - e 4 ) (e 4 - ei)(e 4 - e 2 )(e 4 - e 3 )' 


(B7) 


(B 8 ) 


The numerical integration over momenta has to be 
done with care due to singular denominators. We found 
that numerical integration performs the best using the 
above forms of AF^ after introducing regularization 

Re ( e 3 -eb»r) ’ with r = 10_15 usin § Matlab 
3d integration routine. For details and checks see [28] . 

In 3D, there is a possibility of a non-coplanar inter¬ 
action diagrams. They obtain if there are non-trivial 
{qi,..., q 4 } that add up to 0. Such diagrams exist for ex¬ 
ample for FCC lattice, which has reciprocal BCC. There 
are 8 BCC reciprocal vectors, which can be split into two 
distinct quadruplets (tetrahedra). Each has total 4! = 24 
multiplicity. In the case of FCC, by symmetry all dia¬ 
gram have the same value and each has the same look 
identical to D above. Lets name it D\j. Then, in Eq. 
@> w = 24 D\j. The relative magnitude of coplanar and 
non-coplanar terms is obviously important. In Figure [ 2 ] 
we plotted u(a)/S and in Figure [3] — w/ 24. 

Appendix C: Variational crystalline states 

We have considered the following variational states (N 
is the number of ±q^ pairs, see Figure [5J: 


• Columnar: N = 2. 1 neighbor at optimal anlge 

&min 

• Rhombohedral: N = 3. 2 neighbors at optimal 

angle C^rnin 

• BCC lattice (FCC reciprocal): N = 6. 4 neighbors 
with a = 7r/3, 1 with a = 7r/2; 

• FCC lattice (BCC reciprocal): N = 4. 3 neighbors 
with a = cos -1 (1/3) 

• iQC: N = 6. 5 neighbors with ~ 63.4° 

• Edge-icosahedral (momenta are the edges of icosa¬ 
hedron - favored by cubic interaction which we ne¬ 
glect): N = 15. 4 neighbors with a = 60°, 4 neigh¬ 
bors with a = 72°, 4 neighbors with a = 36°, and 

2 neighbors with a = 90°. In the energy there are 
non-coplanar terms present; we did not include this 
contribution since the energy of this states is rela¬ 
tively too high (due to many suboptimal angles a) 
and the non-coplanar contribution is weighted by 
small factor N~ 2 . 

• Decagonal (same as iQC, but with one vector pair 
missing): N = 5. 4 neighbors at icosahedral angles 


• Smectic or stripe: N = 1 . 
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FCC 






iQC 


Edge-icosa 


decagonal 


dodecahedral 


FIG. 5. Primary reciprocal lattice shells of variational states that we considered. Reciprocal vectors are indicated by red lines 
that start from red sphere and end at golden spheres. Within assumptions of Weak Crystallization theory, only such vectors, 
which all have equal length, contribute to crystallization energy. Clarifications: (1) vertices of perfect octahedron represent the 
primary reciprocal shell of a Cubic state; general rhombohedral case can be obtained by linear dilatation or contraction along 
the axis connecting centers of opposite triangular faces (black line). (2) Reciprocal vectors of the edge-icosahedral state are, as 
the name implies, edges of icosahedron, shown in red. 


• Dodecahedral in momentum space: TV = 10. 3 

neighbors with aq ~ 41.8°, 6 neighbors with ~ 
70.5°. 

• Hexagonal: N = 3. 2 neighbors at 60°. 


Appendix D: Splitting peaks 

Here we show that splitting of one Bragg peak into a 
pair is unfavorable. This is an immediate consequence of 
u(a) being smooth as a —> 0, as is the case for electron- 
mediated and local interactions. Indeed, assume that 
there is an energetically favorable (possibly multi-q) con¬ 
figurations with a spot at qo with amplitude p qo . Now, 
suppose we split it into two at q^ and qg, both approx¬ 
imately equal to qo- To keep the interaction with the 


other momentum components unchanged (we assumed 
it to be optimal), we need || 2 + |p q "| 2 = \Pq 0 \ 2 - 
That keeps the second order (r) and the interaction 
with distant q components intact. However, instead 
of the original self-interaction we now have uo\pq 0 \ 4 

u 0 \ P< \ 4 + 4«(a)|p q , | 2 |p q ,,| 2 « W0 (M 2 + K'l 2 ) 2 + 

2u(a)\p CL ' o | 2 |p q " | 2 - Hence, the energy goes up, and split¬ 
ting is not favored for uq > 0. Indeed, the crystallization 
simulations starting from random initial conditions show 
the extinction behavior: large Bragg peak suppresses its 
smaller neighbors, leaving in the end only a small number 
of spots that correspond to a (q)crystal. 
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Appendix E: Distorted iQC state 

To explore the stability of iQC state with respect to 
distortions away from perfect icosahedron, let us expand 
the interaction energy in the vicinity of the iQC: 

Eint = u { a ij)\Pi I \Pj I • 

i<j 

For the sake of argument will neglect the fact that 
the amplitudes of the order parameter can also react to 
distortions - this will only further lower the energy of 
distorted state. Then, defining Sij = ctij — <ao, 

Eint = Eq + rd(o'o) ^ + 0.5rt // (ofo) ^ ^ Sfj + ... 

i<j i<j 

Now we can define convenient coordinates for the 
Bragg peaks on the sphere, and explore whether the 
energy can be lowered by a distortion. Both the first 


and second derivate terms define quadratic forms with 
non-negative eigenvalues (due to the nonlinear depen¬ 
dence of Sij on local coordinates, even the first order 
term produces quadratic form upon expansion). Out 
of 12 total eigenvalues, the quadratic form of YU<j $ij 
has only 4 non-zeros; in contrast YU K j $ij ^ as on ^y 3 
zero modes that correspond to rigid global rotations. 
When put together, for u'(a o) < — (2/3)u"(ao) nega¬ 
tive stiffness modes emerge, signifying distortive insta¬ 
bility of icosahedron. The strongest instability occurs 
at the largest possible quasimomenta ±47r/5 (see Math- 
ematica code attached to EH). At the critical point 
u'(a o) = — (2/3)iz"(ao)? f° ur zero modes simultaneously 
appear, forming a flat zero-frequency band as a function 
of quasimomentum on icosahedron. 

Hence, the conclusion is that even if u(a) reaches the 
minimum at non-icosahedral anlge, the iQC remains (at 
least) locally stable for u' > 0 (“tensile strain” between 
Bragg peaks), and even for “compressive strain” it re¬ 
main stable until a critical value of negative u' is reached. 


